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The Balitsky-Kovchegov QCD equation for rapidity evolution describing saturation effects at high 
energy admits universal asymptotic traveling-wave solutions when the nonlinear damping becomes 
effective. The asymptotic solutions fall in universality classes depending only on some specific 
properties of the solution of the associated linear equation. We derive these solutions for the recent 
QCD formulations of the Balitsky-Kovchegov equation with running coupling constant obtained 
from quark-loop calculation. While the associated linear solutions depend in different ways with 
observables and higher-order effects, we show that the asymptotic traveling-wave solutions all belong 
to the same universality class whose solutions are given. Hence the influence of saturation stabilizes 
the QCD evolution with respect to higher order effects and leads to universal features at high enough 
rapidity, such as the form of the traveling waves, the intercept of the saturation scale and geometric 
scaling in square-root of the rapidity. 



I. INTRODUCTION 

In the leading- logarithm (LL) approximation of high-energy (high-density) QCD, the evolution with rapidity Y 
of deep-inelastic scattering cross-sections is driven by the nonlinear Balitsky-Kovchegov (BK) equation P, [3, 0). 
This equation is supposed to capture essential features of saturation effects in the "mean-field" approximation where 
fluctuations (or "Pomeron-loop" effects Q) can be neglected. More specifically, for the dipole-target amplitude in 
2-dimensional transverse position space jV(x, y, Y), it reads 

cW(x, y, Y) = £ J d 2 z CT^p {A^(x, z, Y) + Af(z, y, Y) - AA(x, y, Y) - AA(x, z, Y) x Af(z, y,Y)} , (1) 

where \xy\ 2 = (x— y) 2 , x and y are the transverse space positions of the quark and antiquark constituting the QCD 
dipolc. In Eq.|T|), the coupling a = a s N c /i: is kept fixed, since the derivation of the equation is made at leading 
logaritmic approximation, and thus zeroeth perturbative order in the coupling. Indeed, this nonlinear equation 
corresponds to resumming QCD "fan diagrams" in the LL approximation [IJH. 

In the 1-dimensional case, without impact-parameter dependence (e.g. for large targets), 7V(x, y) = Af(\xy\), the 
equation acquires a simple form in the Fourier-transformed momentum space, namely |3J 

a^d Y N(L, Y) = X (-d L ) N(L, Y) - N 2 {L, Y), (2) 

where L(k) — log(fc 2 /A 2 ), and A is here an arbitrary scale. The Fourier-transformed amplitude 

N(L,Y) = f° fei J {k\xy\)M{\xy\,Y) (3) 
Jo \ X V\ 

can be related by convolution (see e.g. [5() to the unintegrated gluon distribution in the 1-dimensional space in 
transverse momentum k = |k|. In this LL approximation the characteristic function of the kernel has the standard 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) form (fj, [7|,l8|, namely 

X (7) = 2^(l)-#y)-^(l-7) ■ (4) 

Quite recently, the extension of the BK equation to running coupling has been the subject of interesting theoretical 
studies. More generally, the extension of the saturation formalism beyond LL order is of present concern. 

The goal of the present paper is to apply the traveling-wave method to the problem of solving the nonlinear BK 
equation with running coupling. This method allows to find asymptotic solutions to nonlinear equations of the BK 
type and to discuss their universality properties. 
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Indeed, we will consider the recent theoretical advances towards a full QCD formulation of the BK equation at 
next-to-leading logarithmic (NLL) order in the coupling constant and beyond. This implies in particular taking into 
account in a proper way the running QCD coupling and NLL and higher order contributions to the kernel. The first 
part of our study is devoted to the equations incorporating the terms determined through the quark-loop contributions 
to the BK equation [TJ, [TTJ, [ll[ (see also related older [l2[ and recent [ID, 0, EH references). Then, we consider the 
full NLL contributions to the kernel [l6L [13] and the higher-order contributions implied by the renormalization-group 
constraints [H, H, 0, HI • 

Since a few yearSjthe solutions of saturation with running coupling have been considered in the literature. Already 
in the papers [22|, [23J , the dominant asymptotic term of the solution has been obtained from the knowledge of the 
BFKL equation. Using this equation with absorbing boundary conditions, the subasymptotic term could also be 
derived and extended to NLL BFKL kernels [H|. 

Our tools for the analysis rely on the powerful traveling-wave method. Indeed, it has been remarked that, contrary 
perhaps to naive expectation, the nonlinear character of the BK equation leads to some simplification w.r.t. the linear 
BFKL equations. Interesting universal properties, i.e. asymptotic solutions which do not depend either on initial 
conditions, on the precise form of the nonlinear terms or on details of the linear kernel. This is tightly linked to a 
mathematical property of a well-known class of nonlinear equations [26l (2?J which admits asymptotic solutions in 
terms of traveling-waves |28l. l29l . [3fJ ■ The relevance of this class of equations and traveling- waves for the BK equation 
has been raised in [HI, l32|. 

There exist previous analytical results relevant for our study on the BK equation with running coupling. Heuristi- 
cally, the running coupling has been considered together with the 1-dimensional BK equation @ by substitu ting the 
fixed coupling in front of the equation with a running coupling a(L) depending on the gluon momentum scale |24Ll32^. 
One writes 



i(L) l d Y N{L, Y) = X (-d L ) N(L, Y) - N 2 (L, Y) 



(5) 



We define L(k) — log(fc 2 /A 2 ), where A is now the QCD scale. One considers the LL kernel function as in Eq.(|3]) and 
the running coupling reads 



a{L) = - , 



11N C - 2N f 
l2N r 



(0) 



In this case, the asymptotic universal solution has been found, either by approximation of the linear regime with 
absorbingconditions [2J, [25[ or directly [HI, HH from the traveling wave method of solving equation §5§ ■ 

In Ref . |9| , the quark-loop contributions to the small- a; evolution has been computed for the dipole-target amplitude 
in transverse-position space. Through a separation scheme between the running coupling at next-to-leading order level 
and the kernel at NLL level, an evolution equation has been derived. For the 1-dimensional amplitude Af(\xy\, Y) the 
correspondingly modified BK equation (see Q, formula (52) where the equation is in fact written using traces over 
Wilson-line operators, easily translated in terms of Af) reads 



'a(\xz\) \ If a(\yz\) 



\xy\ 



\xz\ 2 \zy\ 2 



xz 



a(\zy\) 



yz\ 2 \a(\zx\) 



where 



a(\xy\) 



1 



bL(\xy\) 



L(\xy\) = -log(\xy\ 2 A 2 ) 



(7) 



(8) 



Note that the expression ([7]) needs some regularization procedure, since Landau poles appear in the z-integration over 
a(|a;z|) and a(|y^|). We will consider this important point in detail later on for our study. 

In the same time, another consistent derivation of the BK equation with running coupling constant has been 
performed both for the dipole amplitude in position space and for the unintegrated gluon distribution in momentum 
space (icllllj|. The same quark-loop contribution as in [j| has been found, but using a different separation scheme. The 
quark-loop contributions have been splitted differently into a contribution to the running coupling and a contribution 
to the kernel. The equation found in position space is thus different from ((T)). The involved running couplings appear 
to form a "triumvirate" structure in terms of three different scales. More precisely (see [111 ], where formula (48) the 
equation is written with additional factors e~ 5 / 3 absorbed here in A) one writes 



Y \ = [is ^(g 2 ) «(M) 2 ) 

dY {, ) J 2tt a(fc 2 ) 



- 4>(q,Y) + ± 0(|kHq|,y)- 



k 2 



cp(k,Y) 



nonlinear terms (9) 



with 



<f>(k,Y) = a{k 2 ) <f>(k,Y) , (10) 

4>(k,Y) being the unintegrated gluon distribution in tranverse momentum space. The nonlinear terms in Eq.([9]) are 
not explicitely written, but as already known from other studies of the BK equation in momentum space [33| , their 
precise form is not needed for the derivation of the traveling-wave solutions in the transition region towards saturation, 
provided they ensure the necessary damping. 

This "triumvirate" of coupling constants has been advocated long ago from bootstrap properties [12j, while the 



definition (flT)]) has been introduced more recently . Quite interestingly, the consistency of the equation for the 
dipolc amplitude ([7]) with the one for the unintegrated gluon distribution has been established in [ll[. However, the 
corresponding Pomeron intercept driving the leading rapidity behaviour of the linear equation is different in both 
cases at NLL order, by contrast with the identical rapidity behaviour at LL level. The question of the universality 
properties of the asymptotic solutions of these equations (|5I7I9[) when taking into account nonlinear damping terms is 
thus interesting to investigate. Indeed, the dependence on the separation schemes and the dependence of the Pomeron 
intercept on the observable valid for the linear regime at NLL order and beyond have to be reexamined when including 
saturation effects. 

Concerning the scheme dependence related to the rcnormalisation-group constraints, one starts with the expectation 
that the linear regime of the BK equation is driven by the full NLL BFKL kernel already available from the calculation 
of (T5 , [13] • It is known that one introduces a "renormalization-group improved" (RG- improved) kernel at all orders 
of perturbation. One has to get rid of unwanted spurious singularities brought together with the NLL calculation. 
This boils down to the existence of different RG-improved schemes which are equivalent at NLL accuracy but with 
different resummations at higher orders (see, e.g., the S3 and S4 schemes [l8[ and the CCS scheme (20j). 

The traveling-wave method for the BK equation with the NLL and higher-order kernels have already been used in 
Refs. [34l [35j . In Ref . [3~I| . the NLL kernels have been considered together with a fixed value of the QCD coupling. In 
this case also the method works for finding suitable asymptotic solutions, and they depend on both the value of the 
coupling and on the NLL framework. In some sense, our study is a continuation of Ref. [35| where the RG-improved 
scheme dependence is studied for the heuristic equation (j5]). Note that using the fixed coupling in the definition of 
the RG-improved scheme leads to a scheme dependence of the traveling waves [35]. We want to reanalyze the non 
linear pro blem when using a running coupling for the scheme definition, as done for the corresponding linear problem 
in 0,123. 

It is important to note that, for the Balitsky scheme, Eq.(0, as well as for the "triumvirate case", Eq.Q, the 
question of spurious singularities and of renormalization-group improved kernels have not yet been studied. Moreover, 
the regularization of Landau poles may introduce some extra subtcltics in this problem. We will show, by a careful 
study of singularities in the NLL kernels, that such collinear singularities do not change the determination of the 
universal terms in the asymptotic expansion but they are responsible of a slowing down of the convergence of the 
asymptotic expansion and probably introduce stronger nonuniversal subasymptotic contributions. Such subasymptotic 
contributions have been discussed e.g. in Ref. [251] . They certainly deserve a separate further study in the framework 
of traveling waves. 

The plan of the paper is as follows. In section [Til we consider the BK equation ((5|) with momentum-space running 
coupling and recall the known results on its traveling- wave solutions [32t [35j , including a new discussion of the effect 
of NLL kernel singularities on the universality properties. In section IIIII we turn to position space and derive the 
corresponding linear -reg ime and traveling- wave solutions for the Balitsky scheme of Q. In section Hvl we consider the 
"triumvirate " case [10 . 1 1 11 ] , and exhibits its solutions in the same way. In section [V] we discuss the RG-improved 
scheme dependence. In section I VII we summarize our results by concluding that all corresponding traveling- wave 
solutions fall into the same universality class and indicate the corresponding asymptotic predictions. 



II. BK EQUATION WITH RUNNING COUPLING CONSTANT IN MOMENTUM-SPACE 

A. Solution of the linear equation 

Let us first remind briefly the method used in Ref. [32I ] to obtain the traveling- wave solutions of the Balitsky equation 
([5]). Following the general method [3^|, we first write the solution to the linearized version of the equation. In the 
saddle-point approximation, it has the form of a double Mellin transform [13] 



N(L, Y) = J J ^- JV ( 7 , w) exp (-<yL + uY + 



(7) , (11) 



with the kernel dependence appearing through the function 

*(7)= fdixil') , (12) 



7 being an unspecified constant. Indeed, using the saddle-point method for the integration over 7 at large enough L, 
one gets the equation 



_ L+ _L a 7 x( 7 ) = -L + i- X (7) = - I I-! 



or equivalently in the form of an operator acting on the amplitude 

bL dyN(L, Y) = X (-d L ) N(L, Y) . (14) 

Hence the solution (fTTj) verifies the restriction to linear terms of © , within the saddle-point approximation. 

As a next step one performs the saddle point integration over ui in the limit of large Y . The saddle point lu s is 
given by 



bY ' 

behaving like oj s ~ Y"~5 . With this form of uj s the gluon density is given by 



X{1) (15) 



N(L, y) ~ / %L AT ( 7 ) exp(- 7 L + fi( 7 )t), (16) 

' 27TI 



where the time variable is interpreted as t = yY and 



o ( -i = A /|x(7). (17) 



B. Traveling- wave solutions 

A critical group velocity (defined as the minimum of the phase velocity in the wave language) is obtained as 

Vg = n{ lc )/ lc = n'( 7c ) . (is) 

However, j c determined in such a way still depends on the arbitrary constant 7. Thus, requiring v g to be independent 
on the choice of 7 means dv g (^)/d^ = 0, which in turn gives dv g ('j c )/d r y c = since the dependence of the velocity on 
7 comes through 7c only. Applying this condition to Eq. (|17p one gets 



dVg( 1 ) d 1 = =0 Vg = A — , (19) 

d~i c y b 7c 

eliminating all dependence on the arbitrary constant 7 in the definition (fT2|) . As a consequence, one finds the equation 
for the critical exponent 

X (7c) 

which is the well-known value (j c — .6275) associated to the saturation solutions using the LL BFKL kernel ((3]). 

Finally, using an ansatz technique borrowed from statistical physics [29j], one finds [32j the result for the gluon 
density 



in agreement with [24]. £1 = 



—2.338 is the first zero of the Airy function Ai(£) and the saturation scale has the form 



up to a non universal multiplicative constant. Hence the two first terms of the expansion a large t of d\og(Q 2 )/dt 
(and thus of the saturation intercept d\og(Ql)/dY) are completely specified. Note that the result (f2~Tj) is known to 
be valid in the transition region towards saturation but not in the deep infrared saturation domain. However, the 
infra-red region is also universally constrained by unitarity, namely N ~ log(Qs/fc). 

The important output of this traveling-wave solution with running coupling constant a(L) is that the amplitude 
(|2"Tj) depends only on the ratio k 2 /Q 2 s {t) at asymptotic values of t. It thus corresponds to a geometric-scaling property 
[3a | of the gluon density as a function of t = VY , instead of Y in the non-running case. 

It is useful to note that the ansatz leads to the solution (|21l22p which do not depend on the saddle-point ap- 
proximation used initially in (jTTJ) . Indeed, all parameters of the ansatz are determined self-consistently from the 
asymptotic analysis of the original nonlinear equation. The saddle-point approximation has been useful to determine 
the appropriate form of the ansatz [32j. 



C. Extension of traveling waves to a cj-dependent kernel 

In the following, and throughout the rest of the paper, we introduce the possibility of dealing with a more general 
kernel x{~&l, dy). This generalization comes naturally from the NLL studies. For instance it has been introduced in 
the w-expansion method of Ref.[19j, in order to define RG-improved kernels. More generally, it is useful in order to 
perform a change of variable ct{L) — > uj in the linearized equation, which is the starting point of the traveling wave 
method. It allows to express the solution of the linear equation in terms of a linear superposition of waves. The 
nonlinear term will then select the critical wave as the traveling wave solution. 

Let us recall for further use the extension [35| of the traveling wave method in the case of a kernel x(— <9l, dy) in 
([5|), namely 

a(L)- 1 dyN(L, Y) = x (-d L ,dy) N(L, Y) - N 2 (L, Y) . (23) 
In double Mellin space, the linear part of the equation reads bLuj — ^(7, a;). Introducing 

*(7,w)= f Vx(Y,w) . (24) 



7 



-V(.I.V) : / |^(7) exp 



(25) 



the saddle-point integration over uj leads to the asymptotic solution in Y 

where uj s is given by the saddle-point equation 

Ybiu 2 ~X{ 1 ,Lu s )+Lu s X( 1 ,u Js )=0 (26) 

and the "dot" means the derivative with respect to uj. 

The solution w s (7, Y) of Eq. (|2"6")) defines the asymptotic properties of the amplitude (f2l>)) . In order to solve it, one 
has to take into account the analytic properties of the kernel integral (j2"4")) , which posesses singularities near 7 = 0, 1. 
The solution will depend on the position of the critical value 7 = 7 C for the saturation solution with respect to the 
singularities. Indeed, as we shall now prove, one has to avoid some well-defined regions around the singularities in 
order to find the universality properties. 

Let us first consider the situation where j c is not in the neighbourhood of the singularities of the kernel. Then, the 
expansion of the integral of the kernel (|24|) and its derivative at uj = is analytic near 7 ~ j c . One can use the Taylor 
expansion around uj = 

X^-±*%&*, (27, 



and by collecting terms with the same powers of u) s , Eq. (|26[) writes 



rfc + ix( 7 ,o) 



- S 2 =^(7,0)-|E^^)T^ (P) (7 ; 0)^| • (28) 



Following [35j , the term in braces will contribute only to subleading nonuniversal terms in the traveling wave expansion 
while the left-hand term with the second derivative w.r.t. u> can be absorbed by a translation in Y. This term (and 
the higher derivative terms) depends thus on initial conditions, and we have the universal asymptotic behaviour at 
large rapidity 



- - w 

Hence the previous relations (|15H20[) remains valid once substituting x{l) ~ ¥ x(7> w = 0). 

Now, we have to take into account the singularity regions, in order to discuss the validity range of the previous 
derivation. The singularities of the function X("f, to) depend crucially of the use of a renormalization-group improved 
kernel or not, due to the presence of spurious singularities in NLL contributions of the kernel, which have to be 
cancelled by the renormalization-group improvement scheme. We shall however consider also the case when spurious 
singularities are still present, by comparison. 

• Renormalization-group improved kernels. The behaviour of the kernels x near the singularities are simple poles 
(apart possibly from the quark contributions, which we shall not consider in this derivation). By integration, see 
(|2"4"1) . this leads to logarithmic singularities in the function Xfr, w) and thus in Eg . (f2"6")l . Hence, the singularity 

region to be avoided is 7c , 1 — 7c < O (u> s log \u> s \) = O (log Y /y/Y^j . Inside this region, the expansion ([28]) and 

thus the universal hierarchy in uj s is spoiled. 

• NLL kernels with spurious singularities. The NLL kernels x contain singularities up to triple poles. By the 
integration in (|24[) . new single and double-pole singularities appear at next-to-leading order. Hence, the singu- 
larity region to be avoided becomes 7 C , 1 — 7 C < O . Inside this region, the expansion |28|) 

and thus the universal hierarchy in tv s is spoiled. Thus the useful region of analyticity requires a much higher 
rapidity to be valid. 



This discussion shows that the universality properties will show up at different scales in rapidity, depending on the 
singularity structure of the NLL kernels. In fact, if one considers NLL contributions keeping spurious singularities, 
subasymptotic nonuniversal terms may be dominant up to high rapidities. Hence, the singularity structure of the 
NLL kernels will affect the subasymptotic behaviour. 

We note that the running coupling case, introducing the integration (|24j) . softens the singularities. It may also hap- 
pen that the w-expansion itself softens the singularities (even the spurious ones, if present), hence possibly improving 
the rapidity range of universality properties. This, and the question of quark contributions deserve a more detailed 
study. In the present work we will focus on the case where the critical j c is well inside the analyticity region. Since 
7 C takes the same universal value 7 C = .6275 as for the LL kernel, the universality properties will appear when the 
region of singularities is avoided and thus for large enough rapidities. The price to pay for the eventual presence of 
spurious singularities is that the convergence to the universal features may be sizeably pushed away in rapidity. 

In order to finally get the universal asymptotic terms of traveling-wave solutions, one has to express the linearized 
version of the BK equation (|23[) retaining only terms in the kernel relevant for the asymptote analysis. By expanding 
around u) = 0, 7 = 7 C up to the second order, one gets 

Y t d t N = {- b ~ v * gdL + 1 - x »(dl+2 lc d L+1 t) 

+ i * * ~ 1 X' fcA - i x' 7c* + ^ * (fl? - i ft) } , (30) 

where the "prime" is the derivative with respect to 7. 

The terms in the first line of p0|) correspond to the same expansion as for the w-independent case and they 
contribute up to order Y" -1 / 6 , which is enough to determine the two first asymptotic terms of the amplitude and 
the saturation intercept. The second line contains new terms, with derivatives in to, which contribute only to higher 
order. 



It is interesting to note however that the first dominant term in the second line of (|30p contributes with order y -1 / 3 . 
The nonuniversal contributions are expected to start at order y -1 / 2 , since they correspond to a shift t — * t + to which 
results in nonuniversal terms of order y -1 / 2 . Hence the first term in the second line of (|30p gives a contribution 
related to the w-dependence of the kernel, since it depends on \ which could be the remaining track of NLL effects in 
the universal traveling-wave solutions. 

All in all, one finds the same solutions as (121I22P with the substitution x{l) ~ > xil,^ = 0). 



III. BK EQUATION WITH RUNNING COUPLING FOR BALITSKY'S SCHEME 
A. Balitsky's formalism and regularization 



Let us now derive the traveling- wave solutions for Eq.©. Following the abovementionned method, one first look 
for the solution of its linear part. However, in order to give mathematical and physical consistency, one has to 
introduce some regularization procedure of the integration over z in position space in order to avoid the Landau poles 
at L(\xz\) = and L(\yz\) = 0. 

Starting with solving the linear part of ([7]) 



bL JLtf(\ xy \,Y) = I ^[M{\xz\,Y)+N{\zy\,Y)-M{\xy\,Y)\ 



dY 



R 



2tt 



\xy\ 



1 log(M 2 /M 2 ) 1 log (\xz\ 2 /\zy\ 2 ) 



\xz\ 2 \zy\ 2 



\xz\ 2 log(|xz| 2 A 2 ) \zy\ 2 log(|zy| 2 A 2 ) 



(31) 



where the subscript 1Z denotes the regularization procedure, we insert the Mellin transforms M (7, Y) of the dipole 
amplitude in position space. 



Af(r, Y) = 



2iri 



dr, 



2m I 2ni 



(32) 



Applying the Mellin-transform ([32| to the right-hand side of (|3"T|) . one obtains the diagonalized form of the linear 
kernel 



*r(7,a(M)) 



d 2 ^ 



\xz\ 

\xy\ 



2o 



mi 

\xy\ 



2"; 



^(S) , 1 lo s(&P) 



\xz\ 2 \zy\ 2 \xz\ 2 log(|a;z| 2 A 2 ) \zy\ 2 log(|zy| 2 A 2 ) 



The first term in braces corresponds to the LL kernel The two other terms give an additional contribution to 
higher orders of perturbation which by symmetry reads 



d 2 ^ 

iZ 27T 



\xz\ 

\xy\ 



2~ 



+ 



\m 

\xy\ 



2j 



1 log(|^| 2 /|zz| 2 
\xz\ 2 log(|a;z| 2 A 2 ) 



(34) 



In order to regularize (|34p . we choose to introduce a truncation of the perturbative expansion as a function of 
a(|xj/|) = [— 6 log(|xy | 2 A 2 )] - 1 . We will discuss later the regularization dependence. 

Using a rescaling variable A = (x — z)/\xy\ and going to the complex A-plane one writes 

j / (1-A)(1-A) 

^ [W + ((1-A) (1 -A)r-1] 



' /WA [(AAf + ((1-A)(l-A)f-1] ' 



A) 



1Z 2in 

dXdX 
2in 



XX log(AA) - [ba{\xy\)]-i 



[(AA)V((1-A)(1~A)) 7 -1] ilog 



AA 



(1-A)(1-A) 



A' 



J2(ba(\xy\)) n \og(XX)Y 

n=l 



(35) 



where N defines the level of truncation. For instance, N — 1 corresponds to the contribution of the running coupling 
to the NLL term of the kernel. Note that the only remaining scale dependence comes from a(|a;y|). 



It is then possible to express analytically the expansion (|35p term-by-term. One writes 

X^ a \l,a(\xy\)) = X (7) +<*(\xy\) Jm^ (-&a(M)^V ' (-^ - jL) h(>y,e,S) , (36) 

n— 1 ^ / \ / 

where the integral -Zi(7, e, S) is given in Appendix El 

In particular, the contribution (N = 1) to the NLL kernel reads 

X ^(7,«(M)) = x(7)+MI^I){^(x(7) 2 + *'(7)-*'(l-7))-^p 

= X(7) + i (X(7) 2 + *'(7) - ~ 7)) - ^p} , (37) 

which agrees with the result of , confirmed by [ll| (up to the inclusion of the factor | in a redefinition of A and a 
change 7 — > 1 — 7 w.r.t. [HI). 

Finally, the saddle-point ([7]) gives, in double-Mellin in position space at NLL order 

u J = -L X (l) + T^X {1) (l) , (38) 
bL bL z 

where 

X W( 7 ) = \ (x(7) 2 + *'(7) ~ ~ 7)) - — ■ (39) 
2 7 

As well-known, the NLL contribution (|3"5)) to the kernel has a spurious singularity in I/7 2 . As we shall now see, 
this will not influence the determination of the universal terms in the traveling-wave solution. However, as explained 
in the introduction, it may delay the convergence of the asymptotic expansion. 



B. Derivation of the traveling- wave solution in position space 

Keeping in a first stage the truncation at NLL level, namely N = 1 in (|35I36[) . the nonlinear equation in position 
space can be formally written 

1- dyAf (L, Y) = X f al (~d L , a{L)) Af(L, Y) - A^ 2 (L, Y) , (40) 



a(L) 



where L = — log[(|xy|A) 2 )]. One should care about the exact meaning of Eq. (|4"0"]) . which remains at this stage formal 
in two aspects. The nonlinear term A/"® 2 coming from the quadratic part of the integrant in pip is not explicitely 
given and the operator \i al {— d^, a(L)) is to be well-defined, since it contains two non-comutative variables. 

Concerning the nonlinear terms, we know from general traveling-wave properties that its precise analytic form is 
irrelevant for the asymptotic solutions. The same remains true with the yet unknown NLL contributions to the LL 
nonlinear term. 

For the action of the kernel in Eq. ([4TJ]) we consider the kernel obtained in diagonalized form (|38l39f> , namely 

X f Q '(7,«(i))=x(7) + 6«(i) X W (7) • (41) 
By solving (|38p as a function of w, one obtains a new w-dependent kernel ^1(7, u>) which leads to the equation 



bL, = « 1(7>W ) = x(7) + ^ (7) 2 2+4 ^^M , (42) 

where we have only kept the physical root. Then Eq. (|40p takes the operator form 

bld Y M(L,Y) = Kl {-d L ,d Y )JV(L,Y) -N® 2 (L,Y) . (43) 

The derivation of the traveling- wave solutions comes from a simple extension of the arguments of subsection III CI on 
the extension of the traveling- wave method to a w-dependent kernel. Indeed, Eq. (|4"3"J) is similar to ([23]) in section [Til 
substituting L by L, see ©. 



Taking into account the same same substitution L — > L and considering formula (|42p. the traveling- wave solution 
is driven by the value of the w-dependent kernel ^1(7,0;) at u> = 0, namely 

«i(7,w = 0) = x(7) • ( 44 ) 

We see from Eq. (|44[) . that for lo small enough, that is for the universal terms dominant at large enough rapidity, the 
NLL contribution will not modify the saturation saddle-point 7 C . 

From this general argument, we see that the traveling- wave solution is given from the LL kernel with running 
coupling constant in position space. By the momentum-position substitution property and the independence of 
the result on the precise form the nonlinear damping (despite the fact that they possess a different form in both 
representations) one is led to the universal form of the solution 



Af(L,Y) = const -Y^-Ai |-21og [\xy\Q s (VY)] y^^ j + & j ' (\^y\Qs(VY)) **' , (45) 

where £i,7 c are the same as in formula (f2"Tj) (note that (jlS)) , as well as (|21[) . is not valid in the infra-red. However, 
the infra-red region is also universally constrained by unitarity in position space, namely Af ~ 1). The saturation 
scale is then identical to (|22| up to a non universal multiplicative constant, meaning that the saturation intercept is 
predicted to be the same as (|2"21 at large rapidity. 



C. Truncation and regularization independence 

Let us sketch the general argument leading to the independence of the traveling-wave solution w.r.t. the truncation 
and the regularization scheme. For this sake, we shall use a method of "w-expansion" introduced in the derivation of 
the linear BFKL equations at NLL order with renormalization-group constraints [iji HH • Briefly, this amounts to 
transform the perturbative expansion of the kernel into an equivalent expansion in lo. 

The asymptotic universality property can be extended to a given finite value of N in a more general truncation of 
the expansions (|35l36p . Indeed one can write the corresponding saddle-point equation as 

o J = a(L) X B N a \ 1 ,a{L)) , (46) 

where xfi^ili a {L)) is a polynomial of degree N in a(L). This can be considered as as an implicit equation for L as 
a function of lo namely 

L= ^- k n (j,uj) , (47) 

where the function kn verifies thus 

KN (j, u>) = xT {7, oj/kn^, lo)} . (48) 

Then, coming back to the truncation equation (|3"6"|) . one finds that XN al (li a. = 0) = xil)- Then the implicit equation 
(|4"8f leads to kn(j,lo = 0) = x{l)- O ne again falls in the same universality class. Indeed, this is expected from the 
remark that lo and a(L) are of the same order and jointly go to zero at the limit. 

The same result should apply for any consistent scheme for the regularization of the Landau poles, e.g. by freezing 
of the coupling constants beyond some scale \xz\, \yz\ > ?*o in the integral (|31|) . The perturbative consistency of the 
regularization implies the existence of a convergence domain in 6l,lo near zero. In physical terms, the traveling-wave 
solutions are valid in the domain of high enough rapidity to have a perturbative saturation scale and driven by the 
ultra-violet behaviour. Hence the nonperturbative effects related to the regularization should not be relevant. This is 
a natural constraint on the regularization scheme. 

The regularization, which corresponds to a given kernel Xtz* 1 w ^ then give rise to an implicit equation 

kr(7, w) = Xk 1 {7, uj)} , (49) 



defining a new w-dependent kernel k-r. By the same arguments, one again expects ^7^(7,0; = 0) = x{l)- 



D. Traveling- wave solution in momentum space 



Now we shall demonstrate that the traveling-wave solution obtained in position space in the preceeding subsection 
is in fact in the same universality class as the solution found for the initial BK equation (J5j) with running coupling 
in momentum space. Since a Fourier transform of the asymptotic position-space solution (|21|) is not mathematically 
justified, it is convenient to start directly by transforming the equations from position to momentum space and then 
find the asymptotic solutions. 

Due to the previous universality properties in this space, the problem reduces to the solution of the initial equation 
([7]) when keeping only the LL kernel: 



bi JLm{\x V \, y) 



d 2 z \xy\ 2 
2ir \xz\ 2 \zy\ 



[Af(\xz\, Y) + Af(\zy\,Y)-Af(\xy\,Y)-Af(\xz\,Y) Af(\zy\,Y)] 



(50) 



where the coupling depends only on the parent-dipole size \xy\. 

Then, performing the Fourier transform (J3]) on the two sides of the equation gives 



I" j-y h(Mxy\) b~L -jLtf<\ xy \,Y) = x(-d L )N(k, Y) - N 2 (k, Y) . 
Jo \ X V\ 01 



Using the identity 



one arrives at the equation: 
bL^N(k,Y)- 



d\xy\ 

\%y\ 



bL(\xy\)=bL(k)-b\og(\xy\ 2 k 2 ) 



d 



3 (kr)b\og(\xy\ 2 k 2 )—Af(\xy\,Y) = X {-d L )N{k, Y) - N 2 (k, Y) , 



(51) 



(52) 



(53) 



where an extra term appears with respect to (0, depending still on the derivative of the position-space dipole 
amplitude Af. 

Denoting N the Mellin transform of AT(\xy\, Y) and N, the Mellin transform of N(k, Y), they can be simply related 
(see appendix IB"]) by 



iv( 7 ,y) = 2 2 ^ 1 f ^^^( 7 ,y). 

For the derivative term, by insertion of Mellin-transform, one writes the relations 

«JoW^)lo g (|^)A|^ e -,^ (7iy) 



(54) 



= b J |L ^tf(i,Y)A^ J" d\xy\Mk\xy\)togQxy\ 3 h?) {xy^ 1 
dry d 



2ni dY 



Af( 7l Y) 



A 



2- 



2 2 7 - 



duJo(u)u 27 1 = 

1 r( 7 ) 



r(i- 7 ). 

Using now the relation (|54p between Mellin transforms, one finally finds the simple expression 

^e-M7)^V(7,r), 



with 99(7) = ^ log 



927-1 r( 7 ) 
z r(i- 7 ) 



#(7) + *(1 - 7) + 2 log 2. Hence, ||53J) becomes 



bL—N{k, Y) = X (-d L )N(k, Y) + bv(-d L )—N{k, Y) - N 2 (k, Y) . 



Following the same derivation as previously, the linear solution of the linear regime writes 

^ (L ' y) = J Si J ^^o(7^) exp (-^L + u;Y + ±-X( 7 )+J\ 1 , V ,( 7 '] 



(55) 



(56) 



(57) 



(58) 



where A (7) is the function considered in the LL case, see (|12p . This solution can thus be identified as a change 
of the kernel eigenvalue x(l) ~> xili 1 ^) = xil) + buxp^). Hence, the universality class of the non-linear equation 
(|57| is determined by xili^ = Q) = x(l)- This result can also be obtained by noting that the term J7 1^7' 93(7') = 

log 2 27-2 ^ p[^jp|^I^ can be absorbed in a redefinition of the impact factor Ao(7,w). As easy to realize from the 

initial saddle point equation (|13[) , a contribution of the impact factor may only consist in a change of the non universal 
reference scale Qq. 

The final conclusion is thus that the traveling- wave solutions of the Balitsky equation ([7]) are in the same universality 
class as the original BK equation ([5]) with LL kernel both in position and in momentum space. Note that this result 
extends in a non trivial way to running coupling a property which was easy to derive when the coupling is fixed and 
thus factorizcd from the Fourier transform. 



IV. TRIUMVIRATE OF RUNNING COUPLINGS 



Let us now apply the traveling- wave method to the evolution equation ^ . As previously, one has to consider first 
its linear part. As in section ITTT1 one has to introduce a regularization procedure to avoid the Landau poles, except 
that the integration is now done in momentum space. Let us choose the same type of truncation 1Z as for the previous 
equation: we expand the Landau poles denominators and truncate the resulting series. 

Starting with the linear part of Eq.© 



9 fck,Y) = 



dY 



d 2 q a(<? 2 )a((k-q) 2 ) 



2n 



a{k 2 ) 



0(k,Y) 



we rewrite the triumvirate of running coupling 

a{k 2 ) 



a(q 2 ) a((k-q) 2 ) _ ^ log 2 (^) 



a(k 2 



a(k 2 



log(£) log(^) 



l + fe(fc 2 )log(|r) l + 6a(fc 2 )log 



-(k-q) 



(59) 



(60) 



Inserting ([60]) and the Mellin representation <fi(k, Y) = J 
the diagonalized form of the regularized kernel Xti" 1 



2m 



0(7, Y) in the right-hand side of (|59|) . one obtains 



xk r \iMk 2 ))= ^ 



CTq I (k-q) 2 \k 2 J + q 2 \ fc 2 J 



k 

9 2 (k-q) 2 



In 2n |(l + &a(fc 2 )log(^)) (l + 6a(fc 2 )log(^^)) 
Using a rescaling A = q/fc and going to the complex A plane, one writes 



(61) 



dXdX 



n Am [l + 6a(fc 2 )log(AA)] [l + ba(k 2 ) log((l — A)(l — A))] 



(AA) 
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(1-A)(1-A) (AA) 
E (^(fc 2 )) n+m / ^— log"(AA) log m ((l - A)(l - A)) 



N N — r 
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(AA) 
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(1-A)(i-A) 



((1-A)(1-A))- 
(AA) 



(AA(1-A)(1-A)) 



(62) 



where N defines the level of truncation. For instance, N = corresponds to the LL BFKL kernel, and N = 1 to the 
NLL term of the kernel. It is then possible to express the expansion (|62p analytically term-by-term. One formally 
writes 



N N- 



QTl QT, 



n— m— 



(63) 



where the generating functional 12(7, e, S) is given in Appendix [Cl 
At the NLL level one finds 



y(V) 2 3 



o 



(64) 



Once going from the function <fi (see (jlOp ) to the unintegrated gluon distribution 
factor 3/2 into 1/2, this result gives back the one found in 
Eq.Q then reads in double Mellin space at NLL order 



1 / x 1 1 
bL x[l) + bL 2 X 



(i) 



(7) 



which results in changing the 



(65) 



where 



X 



th (i) (7) = Xizt + 3 ($/( 7 ) _ $'(1 _ 7)) 



2 2 

We note that this result is different than for the dipole amplitude, (|39p. as discussed in 



(66) 



Keeping in a first stage the truncation at NLL level, namely N = 1 in 
space can be formally written 



the nonlinear equation in momentum 



1 



a(L) 



Y N(L, Y) = X l"(-d L ,a(L))N(L, Y) - N® 2 (L, Y) 



where 



X J"(-d L ,a(L)) = X (-d L )+ba(L) X q 



(i) 



(-Ol) 



(67) 



(68) 



with the nonlinear contribution in momentum space denoted by N® 2 (L, Y). Again and also in previous QCD traveling 
wave studies at nonzero transverse momentum [33j the precise form of N® 2 will not matter for the universal behaviour 
of asymptotic solutions and thus need not be explicitely derived. 

As in the previous section, one has to solve (|65p as a second-order equation for L. One finds 



bLu 



X ( 7 ) + v/ X ( 7 ) 2 + 4^ X T " «(7) 



(69) 



As explained in the introduction, and valid already for the Balitsky scheme, the spurious collinear pole contained 
in the NLL term of (|68p will not change the universal terms, while it may delay the convergence of the asymptotic 
expansion. 

The whole derivation of traveling wave solutions parallels the one described in the previous section, except for the 
transformation from position to momentum space. Hence finally the traveling wave solutions are directly expressed 
as in (|21I22|) . 

Once again, the generalisation to a fixed perturbative truncation is expected to hold. The universality class of a 
regularized BK equation with the triumvirate of running couplings is the same as the equation including the running 
coupling with the "external" gluon momentum (J5]). Thus, including the effects of the transition to saturation leads 
to an unification of the asymptotic solutions, irrespective of the differences at the level of the solutions of the linear 
regime. 



V. RENORMALIZATION-GROUP IMPROVED SCHEME DEPENDENCE 



Let us recall why the solutions of the BFKL equation with NLL corrections require some special treatment and why 
this introduces a NLL scheme dependence. Indeed, although the NLL corrections are known [l6|, [lTj they turn out to 
be negative and large. This is due to spurious singularities which are in contradiction with constraints coming from 
the QCD renormalization-group properties. Indeed the cancellation of these singularities can be obtained by suitable 
contributions at higher orders. However, since these contributions are not yet calculated, there is a need for convenient 
schemes ensuring the compatibility with the renormalization group, which are not uniquely defined. If one admits 
that the linear part of the various nonlinear evolution equations we are studying will be driven by the corresponding 
BFKL kernels, we have to take into account this scheme dependence in our discussion of universality properties. A 
first discussion of the scheme dependence can be found in ref.[35j], where different schemes were discussed for the 



kernel of equation © . The aim of the present section is to extend the discussion to the equations considered in the 
previous sections and to consider the running coupling also in the scheme definition. Note that the RG-improvement 
plays already an important role in the convergence aspects of the universality properties, see section II-C. 

For our discussion, one could consider at least two different classes of schemes, following Ref. [H[. In a first class, 
containing e.g. the schemes of Refs. [2(J HH, the dependence of the scheme on the higher orders of the perturbation 
expansion has been expressed through the ^-dependence, following the method of [l9(. In this case, the higher order 
resummation appears in the kernel only through the dependence over the two Mellin variables, i.e. x{l — lu = 3y), 
where lo = 0(6l) drives the higher-order corrections. For this class of two- variable NLL kernels, the general argument 
that the traveling- wave solutions are driven by x{li U} = 0) = x(t) applies, see section II-C. 

In a second class of schemes such as S3 and S4 [l8j], the analysis of [35[ leads to results depending on the value 
chosen for the fixed coupling appearing in the definition of the renormalization-group scheme. Indeed, the traveling- 
wave solutions are driven by the kernel x(7j^ — 0,a 7^ 0) ^ x{l)- Hence, even if the form of the solutions is similar 
[HI, the critical parameters are different and thus the traveling- wave solutions appear to be scheme-dependent. One 
practical question, for instance, is to know what are the predictions for the traveling-wave solutions when choosing 
the coupling defining the scheme varying with the momentum scale, see e.g. [38j . Indeed, one would prefer a unified 
treatment of the running coupling effects taking into account the running coupling in the RG improved kernel itself. 
We shall show, that in this case, and on variance with the results at fixed coupling for the regularization, the same 
universality class defined by the LL kernel (with running coupling) is recovered. 

For this sake, we can consider in the same way all previous equations, either that corresponding to the BK equation 
in transverse momentum, (TIB")) , or in transverse position space (f3"5]) . or for the case of the triumvirate ([55]) . using any 
of the renormalization-group improved NLL kernels. Then the corresponding saddle-point equations become 

lu = a(L or L) x NLL {l, u, a{L or L)) . (70) 

Now, in parallel with the previous discussion, Eq. (|70p can be interpreted in all cases as an implicit equation for L or 
L namely 

LorL = ^K NLL (rf,w), (71) 
where the function k nll verifies the implicit equation 

= ^{^-^-^. (72) 

This boils down to reformulate schemes, starting from S3 and S4, by using an appropriate w-expansion (l9| . These 
schemes are now in the same first class discussed above. 

Then, the condition k nll (^, uj — 0) = ^(7) is expected to be verified in all cases by perturbative consistency, and 
thus by the same argument, one again falls in the same universality class. Indeed, if we use the w-expansion method 
to redefine the schemes S3 or S4, it is possible to show that K NLL (-f,u = 0) = xil^ = 0, a = 0) = x{l)- All other 
schemes verifying this relation will again fall in the same universality class, defined by the traveling-wave solutions 
(|21|22p . The question remains for further study whether other types of NLL schemes, such as the one of Ref.[H| can 
be given a similar treatment. 



VI. SUMMARY AND PREDICTIONS 



To summarize, we considered three versions of the QCD evolution equation in rapidity in the mean-field approx- 
imation, i. e. the Balitsky-Kovchegov [J @, Q equation ^ , extended to take into account the running coupling of 
QCD. The first one considers the 1-dimensional BK equation for the dipole amplitude in momentum space, with the 
substitution of the fixed coupling by a running coupling in terms of the gluon transverse momentum, see ([5]). The 
two other forms come from recent theoretical calculation of the quark-loop contributions. They differ by a separation 
scheme for higher orders between the running coupling and the kernel. The second equation we consider is written 
for the dipole amplitude in transverse position space , see ([7]) [§]. The last one can also be written for the dipole 
amplitude in transverse position space, leading to a "triumvirate" of running couplings, and differs from the previous 
one. We study it in its equivalent formulation as the evolution equation for the unintegrated gluon distribution in 
momentum space see © [Tol . [ill ] . We have also enlarged the discussion by considering the modified BK equations 
including the renormalization-group improvements of the NLL BFKL kernels [18|, |20j, l2l| . 

Let us summarize our results. The saturation effects which are formulated through the nonlinear damping terms 
in the BK equation leads to asymptotic traveling- wave solutions for all equations. Their behaviour at high rapidity 



is highly universal. More precisely the form of the solution at high rapidity and the two first terms of the rapidity 
expansion of the the saturation intercept d\og(Q s )/dY are the same. Hence they are identical to those found and 
derived in [32l ] for the initial BK equation with running coupling ([5|). 

For establishing this strong universality property, We have examined the different types of higher order dependence 
which modify the solutions of the linear regime 

• Observable dependence: The traveling wave solutions at high rapidity happen to be independent of the considered 
distribution functions, by contrast with the solutions of the linear part of the equations. For instance the dipole 
amplitude and the unintegrated gluon distribution which lead to different rapidity dependence in the linear 
regime for the same scheme have the same saturation intercept and traveling-wave spectrum. This leads to 
an unified asymptotic NLL predictions for observables, independently of its formulation in terms of dipole or 
gluon distributions. 

• Separation scheme dependence: By computing the traveling-wave solutions for both the Balitsky scheme [|| and 
the Kovchegov-Weigert one [l(J EH, we the same results, belonging to the universality class of the simpler 
equation([5]) . 

• Regularization dependence: The definition of equations (|7I9[) requires a regularization, since their consistent 
formulation should avoid Landau-pole singularities. Using a finite truncation of the QCD perturbative expansion, 
we show that the order of truncation, which affects the linear regime, does not change the asymptotic traveling- 
wave solutions. We propose a general argument for general perturbatively consistent regularization schemes. 

• RG-improved scheme dependence: The QCD coupling appears also in the definition of appropriate schemes for 
using the NLL BFKL kernels in order to avoid spurious singularities. The improved NLL kernels with constant 
coupling were shown to give scheme-dependent traveling wave solutions [3J|, even if the overall coupling is 
running [35j . If the scheme is defined using the running coupling through the w-expansion [l9| , we recover the 
strong universality property. 

• Position vs. momentum- space dependence: The universality class is the same when the equations are expressed 
either in position, or momentum-space formalism by Fourier transformation. This shows that the saturation 
intercept is the same. The form of the front is also invariant except for the infra-red (with respect to the 
saturation scale) regions which are not bound by universality properties. We have already mentionned that the 
infra-red region in both representations is also universally constrained, by unitarity, namely M ~ 1 in position 
space and N ~ log(Qs/fe) in momentum. 

• Rapid vs. slow convergence: By discussing the effects of the singularities contained in the NLL kernels, we have 
distinguished between the RG-improved kernels and those containing spurious singularities. Indeed, the rapidity 
region where the universal properties appear is expected to be delayed in the case where such singularities remain. 
It is thus important to use, whenever known, a RG-improved scheme for exhibiting the universality properties. 

From our results, we conclude that the saturation effects on the rapidity evolution of the dipole amplitude or the 
unintegrated gluon distribution with running coupling give a stabilisation of the solutions with respect to higher order 
corrections at high enough rapidity. They appear to be independent of higher order contribution either to the linear 
kernel or to the nonlinear damping terms, however, subasymptotic nonuniversal contributions may be present. 

In mathematical terms one expresses this property as the existence of a large universality class of solutions given 

by mnp . 

At this stage, it is possible to formulate predictions for the behaviour of the exact solutions of the evolution equations 
which can be obtained through numerical simulations. We can list them as follows: 

• The form of the traveling- wave fronts (except in the full saturation region) should converge at high rapidity for 
the solution of the different equations in the same representation space, either in position or momentum. 

• The saturation-scale intercept d\og(Q s )/dY which can be determined from these fronts should converge to the 
same function. 

• Geometric scaling in yY, related to the constant limit of the wave speed should be seen from the solutions, at 
least in the middle of the front. 

These predictions can be tested using numerical simulations of the rapidity evolution corresponding to modified BK 
equations with running coupling. In Ref.[40(, the running coupling was heuristically introduced for models in position 
space and compared with the fixed coupling case. Naturally, they could not include the recent theoretical advances, 
but some of the conclusions are still interesting to quote and seem to fit approximately with our predictions, even if 



the equations are different. The traveling- wave regime is observed at high enough rapidity, the saturation scales for 
running coupling seem to converge, at variance with the fixed coupling case where they depend on the value of the 
coupling. The form of the front verifies geometric scaling. 

While completing this theoretical analysis we were informed of interesting numerical simulations of the relevant 
equations [41| using a regularization scheme based on the freezing of the coupling. The simulations [4l[ seem to 
indicate that the traveling- wave structure is preserved at moderate rapidities and independent of the initial conditions. 
However, the wave parameters seem to differ from the critical ones. We predict that at higher enough rapidity this 
difference will decrease and disappear. Note that arguments based on subasymptotic parametric solutions of the 
nonlinear equations [42j show that the traveling wave structure can be maintained at moderate rapidity but with 
modified speed as observed in [4111. Also the form of the wave front is modified, which could give an explanation of 
the difference already noticed in [40( between the running and non running cases. 

On a phenomenological ground, such a theoretical solution is characterized by a geometric scaling [36| in VY , either 
in momentum or position space. By contrast with fixed coupling which leads to geometric scaling in Y [HI, El, this 
property is characteristic of the BK equation with running coupling [25|, H2] . Geometric scaling in Y or yjY seem to be 
consistent with data on the proton structure functions [37j • However, the precise comparison with the parameters of the 
theoretical traveling-wave solutions show that nonuniversal terms are present, and thus nonasymptotic contributions 
cannot be neglected. Parametric solutions, as in [HJ could help to take into account nonuniversal contributions. 

As an outlook for further directions of theoretical study, we mention the extension of our investigation beyond 
the mean-field approximation leading to the BK equation. For this sake, one has to consider the hierarchy of QCD 
evolution equations including fluctuations, i.e. Pomeron-loop terms, when the coupling constant is running. Also 
the role of impact factors, e.g. for the coupling to the virtual photon in deep inelastic scattering has not been yet 
studied at NLL level. Let us also note that a third (and last) universal term in the asymptotic expansion of the 
saturation intercept, which is known for the fixed coupling case [3~l| . could be the remaining track of the NLL kernels 
in the universal traveling- wave solutions as discussed after equation (|30|) . Finally a specific study of subasymptotic 
terms is deserved. It may lead to traveling-wave properties at more moderate rapidities and thus interesting on a 
phenomenological point of view. 
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APPENDIX A: CALCULATION OF h(j,e,S) 



The integral to be computed in (|35| is 

h( 7 ,e,5) = f^mr +((1-A)(l-A)f -1] 
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Applying the expansion operator (|36|) on the result (|A2j) . one writes 
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At first next-to-leading order, one finds 
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APPENDIX B: MELLIN-TRANSFORMS IN POSITION VS. MOMENTUM 

The Mellin-transforms of the dipole amplitude arc defined in position space as 

and in momentum space as 



' h - ■ fdl f dUJ ^-^Ni^u) 
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Let us relate these Mellin-transformed amplitudes by rewriting the initial dipole amplitude in momentum space ([3]). 
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Hence, by comparing (|B2[) with the last expression in Eq. (|B3p . the two Mellin transformed functions are found related 
by 



N( 1 ,Y) = 2^- 1 T ^- ) M( 1 ,Y) 



(B4) 



APPENDIX C: CALCULATION OF J 2 ( 7 ,e,5) 

The integral to be computed in (f6U)) is 
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The next-to-leading result comes from the expansion 



I 2 (j,e,S) = X {l) + (e + S) 



x^ + |(*'( 7 )_*'(i_ 7 )) 



+ 0(e 2 )+0(ed)+0(S 2 ) 



(CI) 



(C2) 
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